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Reverse Micelles Enable Strong Electrostatic Interactions of Colloidal Particles in 

Nonpolar Solvents 
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While the important role of electrostatic interactions in aqueous colloidal suspensions is widely 
known and reasonably well-understood, the relevance of charge in nonpolar suspensions remains 
mysterious. We demonstrate that particles can have surprisingly strong electrostatic interactions 
in low dielectric constant environments when ions are solubilized by reverse micelles. A simple 
thermodynamic model, relating the structure of the micelles to the equilibrium ionic strength, is in 
good agreement with both conductivity and interaction measurements. Since dissociated ions are 
solubilized by reverse micelles, the entropic incentive to charge a particle surface is qualitatively 
changed, and surface entropy plays an important role. 
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Charge is a powerful tool for manipulating the prop- 
erties of complex fluids. In aqueous colloids, charges 
readily dissociate and play an essential role in determin- 
ing the structure of suspensions and their response to 
electrical, mechanical and optical stimuli. Perhaps most 
importantly, charge stabilizes aqueous colloidal particles 
against flocculation. By contrast, charge is commonly 
thought to play no role in nonpolar environments in ther- 
modynamic equilibrium. However, a close examination of 
existing technologies suggests a richer story. Most stun- 
ningly, by assigning opposite charges to absorbing and 
reflecting particles, nonpolar suspensions have been en- 
gineered to act as electrophoretic ink Furthermore, 
charge control agents have long been added to petroleum 
products to reduce the risk of explosions during process- 
ing and transport 0- There is even some evidence that 
charge contributes to the stability of soot in diesel en- 
gine oil 3]. In these cases, reverse micelles play a cru- 
cial role by reducing the electrostatic penalty of ioniza- 
tion and thereby stabilizing ions in equilibrium 
all cases, it is absolutely essential that ions are made 
soluble in this low dielectric constant environment. Re- 
cently, force measurements have provided the first direct 
evidence of electrostatic repulsions between charged sur- 
faces in a nonpolar solvent . Despite the growing ac- 
knowledgement of the importance of charge in nonpolar 
systems, a useful quantitative framework to describe the 
origin and interactions of charge, comparable to the the- 
ory of Derjaguin-Landau-Verwey and Overbeek (DLVO) 
for aqueous systems [EQ, does not exist. Such an un- 
derstanding is essential to fully exploit the technological 
capabilities of charge in nonpolar solvents and to fully 
comprehend the behavior of these fascinating systems 

In this Letter, we present a physical picture of the ori- 
gins of charge and the interactions of colloidal particles in 
nonpolar solvents with reverse micelles. Our experiments 
show that particles charge and interact strongly with 
screened- Coulomb interactions. Ions reside within the 
polar cores of reverse micelles. Reverse micelles not only 



stabilize charge dissociated from particle surfaces, but 
also charge spontaneously to provide additional screening 
of electrostatic interactions. The confinement of charge 
within reverse micelles qualitatively changes the entropic 
incentive to charging: the configurational entropy of the 
particle surface charge, typically negligible in aqueous 
systems, plays a critical role. 

We explore the interactions of colloidal PMMA 
particles suspended in a mixture of oil and surfac- 
tant. The particles have a radius of 780 nm and 
are sterically-stabilized by a layer of poly(12-hydoxy- 
steric-acid) (PHSA) hairs grafted to their surface. The 
particles are suspended in a mixture of dodecane and 
an ionic surfactant, di-2-ethylhexylsulfosuccinate (AOT). 
Above its critical micellar concentration, AOT forms re- 
verse micelles whose size is very nearly independent of 
concentration, |^, 0|. These micelles are very hygro- 
scopic. To minimize the effects of moisture, we mix pure 
dodecane (99%) and AOT (98%) in a dry glovebox. All 
measurements are either performed in the glovebox or 
with sample cells that have been sealed in the glovebox. 

In pure dodecane, our PMMA particles reversibly ag- 
gregate, as shown in Fig. 1(a). However, when AOT is 
added at concentrations above its critical micellar con- 
centration, the particles disperse. Moreover, particles in 
suspensions with 12 mM AOT rarely approach within 
three particle diameters in each other, as shown in Fig. 
1(b). We quantify this surprisingly strong repulsion by 
collecting a series of statistically independent images of a 
freely-diffusing monolayer of particles confined between 
two microscope slides. The centers of each particle are 
then precisely located to calculate the ensemble- 
averaged two-point density correlation function, g{r), us- 
ing a total of about 40,000 particles. As shown in the 
inset of Fig. 1(c), particles in a 12 mM suspension never 
come within four microns of each other. At larger sepa- 
rations, g{r) rises slowly to a peak near 7 /xm and decays 
to unity at larger separations. In order to minimize the 
effect of the confining surfaces, the surfaces of the sam- 
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FIG. 1: Charge-stabilization of a nonpolar suspension, (a) 
Optical micrographs of PMMA particles in pure dodecane 
and (b) with 12mM AOT, field of view 135 x 108 fim^ . (c) 
Inferred particle pair potential, symbols, and fit to screened- 
Coulomb interaction, line. For comparison, a hard sphere 
interaction is plotted as a dashed line. Inset. Measured and 
simulated equilibrium pair correlation functions, symbols and 
line, respectively. A dashed line indicates contact. 



pie cell are prepared to mimic the surface of the particles 
and interactions are measured in the thickest regions of 
the cell that contain a well-confined monolayer [ijj 

The observed pair-correlation function is consistent 
with the repulsive component of the DLVO theory of 
colloidal interactions. We extract pair potentials by in- 
verting g(r) with liquid structure theory [l^ . Using the 
resulting potentials, 2-D Monte Carlo simulations of our 
experiments generate pair-correlation functions that are 
in good agreement with those measured. A comparison 
of measured and simulated g(r)'s for a sample with 12 
mM AOT is shown in the inset of Fig. 1(c), where data 
and simulation are indistinguishable. The observed po- 
tential, shown in Fig. 1(c), is in excellent agreement with 
the screened- Coulomb form. 
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where C is the apparent potential at the surface of the 
particle. The Debye length, — \/ AirXBriion, char- 
acterizes the range of electrostatic interactions, where 
n.ion is the number density of ions. The Bjerrum length, 
\b — /ATreeoksT, is the separation of two elemen- 
tary charges when their electrostatic interaction energy 
is equal to ksT. A least-squares fit to our data yields 
e|C|/fcsT = 4.4±0.3 and k = (1.07 ±0.05) fim-^. Defy- 



ing conventional wisdom, this zeta potential is as hig h as 
those found in highly-charged aqueous systems |l3l H^ . 
However, the effective charge. 



C(a/AB)(1 + Ka), 
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is 225 ± 15, two orders of magnitude below that measured 
in similarly sized charged particles in water. This charge 
on the surface of particles must be balanced by counte- 
rions in solution to ensure charge neutrality. Using our 
measured value of Z and the number density of particles, 
we calculate the density of counterions to be about 0.5 
/im^"^. However, this ion density is too low to account 
for the observed value of k, which implies a significantly 
larger ion density of 3.3 ± 0.3 ^m~^. Therefore, parti- 
cle surfaces must not be the only source of stable ions in 
solution. 

Indeed, stable ions exist in solutions of AOT in dode- 
cane without particles. While these ions are too small to 
be imaged optically, we can detect their presence through 
their contribution to solution conductivity. Assuming 
that charge is carried by singly charged micelles, then 
the conductivity 
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where an is the hydrodynamic radius of the micelles and 
rj is the solvent viscosity. Our 12 mM AOT solutions 
have a conductivity of 17 pS/cm. Using the known hy- 
drodynamic radius of AOT reverse micelles in dodecane, 
we find an ion density of 2.7 ± 0.4 ^m~^, consistent with 
K — 0.98±.05 /im~^, in good agreement with our interac- 
tion measurements. Therefore, ambient charged micelles 
dominate the ionic strength of the system. However, by 
comparing the number density of ions to the number den- 
sity of micelles, calculated from AOT concentration and 
micellar aggregation number, we find that only one mi- 
celle in 10^ develops a charge. 

To elucidate the origin of this excess charge, we mea- 
sure the scaling of solution conductivity with AOT con- 
centration. Varying the concentration of AOT from 0.8 
to 200 mM, the conductivity increases linearly from 1 
to 500 pS/cm, as shown in Fig. 2. Therefore, the ratio 
of ions to micelles, x = nion / nmiceiie is indepedent of 
the concentration of AOT, and x = 1.2 x 10~^. Only a 
tiny fraction of micelles ionize to form a charge carrier. 
While weak electrolytes, of the form AB ^ yl+ + B~ , 
can exhibit similar ionization fractions, x would scale 
like Instead, a charging mechanism where neu- 

tral micelles reversibly exchange charge though a col- 
lision, A + A ^ A'^ + A~ , leads, through the law of 
mass action, to a ionization fraction independent of the 
concentration of A 0, 0|, as demanded by our data. 
The fraction of charged micelles is determined by the 
difference in free energy between the charged and un- 
charged states, X = 2exp(— AF/2fcBT). If the free- 
energy cost of ionizing a single isolated micelle is gM, 
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FIG. 2: Conductivity of AOT/dodecane solutions without 
particles. Symbols indicate measurements, and dashed curve 
indicates fit to Eq. I^with empirical [AOT]-dependent solu- 
tion viscosity, inset Schematic diagram of two-body process 
that leads to the creation of charge in the bulk of the solution. 



then X — '^Q'y^Y>{~9M /kBT). Comparing this basic ther- 
modynamic result to our measured value of x, we find 
that gM — 12 ksT. 

The measured micelle ionization energy has a simple 
electrostatic interpretation. A singly-charged conducting 
sphere of radius surrounded by a dielectric medium 
has a electrostatic energy Ab/20c Assuming that 

electrostatic contributions dominate the ionization free- 
energy and the core of the micelle is effectively a conduc- 
tor, then our measured value of gM implies a core radius, 
flc — 1.2 nm. This is in good agreement with neutron 
scattering measurements jl7j . Thus, we conclude that 
the equilibrium ionic strength of a solution of reverse 
micelles is determined by a competition between entropy 
and electrostatics. The relative size of the Bjerrum length 
and the micellar core radius determine the charge frac- 
tion. When As ^ flc a small fraction of micelles charge. 
When \b ~ Oc, nearly all of the micelles can charge and 
individual micelles can become multiply-charged. 

Ionic strength has a dramatic effect on the interactions 
between particles. We find that the range of the inter- 
actions is greatly reduced as the concentration of AOT 
increases from 3 to 200 mM, as shown in Fig. 3. At 
the lowest concentration of AOT, particles feel a repul- 
sion greater than fcsT at center-to-center separations less 
than 7 /im or 4.5 particle diameters. At the highest con- 
centration, particles can approach each other within 2.5 
jjLxn.. These data are all fit by the DLVO form with a 
constant apparent surface potential ejCI/fesT = 4.1 ±0.1 
and ranging from 0.2 to 1.4 /im. The ionic strengths 
inferred from fits to interparticle potentials are in good 
agreement with the ionic strengths inferred from conduc- 
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FIG. 3: Variation of interactions with [AOT]. Inferred inter- 
actions for solutions of 3, 12.5, 50 and 200 mM AOT, squares, 
circles, diamonds and triangles, respectively. A hard sphere 
potential is shown as a dashed line for comparison. Screened- 
Coulomb fits to each potential are plotted as solid lines, inset 
Comparison of k as inferred from conductivity and equilib- 
rium structure, crimson squares and black diamonds, respec- 
tively. 

tivities, as shown in the inset to Fig. O However, the 
discrepancy between these two measures widens at the 
lowest concentration of AOT, refiecting the contributions 
of counterions to screening. Interestingly, we consistently 
observe a very small attraction in 200 mM samples. We 
believe that this attraction refiects the underlying struc- 
ture of the reverse micelle solution, which displays a pro- 
nounced increase in viscosity at this concentration. 

Our conductivity and interaction measurements show 
that thermodynamics is an excellent tool for predicting 
the concentration of ions in solution. Thermodynam- 
ics can also describe the charging of colloidal particles. 
In order for a particle to become charged in equilib- 
rium, the energy cost of creating and separating sur- 
face charges and counterions must be overcome by en- 
tropy. This energy cost has three components. As we 
have already shown, the counterion has a solvation en- 
ergy, gM = liksT. Similarly, an ion on the particle's 
surface has a solvation energy, gs- Once the ion and 
counterion have been created, they must be separated at 
an energy cost of e where C is the electrostatic poten- 
tial at the surface of the particle. These costs are over- 
come by gains in the configurational entropy of micelles 
in the bulk and the configurational entropy of charges on 
the particle surface. When the particle charge increases 
by one, a neutral micelle in the bulk is converted to a 
charged micelle. Thus, the change in the free energy is 
kBT\\i.[nion/2nmiceiie] = fcBTln[x/2], where we have as- 
sumed that only a tiny fraction of micelles are charged, 
and that the entropies of the charged and neutral ionic 
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species have the same form as an ideal gas. If a particle 
has N chargeable sites and a charge of Z with Z <^ N, 
then the surface entropy's contribution to the chemical 
potential is kBTln[Z/N]. Thus, the total chemical po- 
tential is 



= 9S + 9M 
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In equilibrium, /i = 0. In the limit where excess ions 
dominate over counterions, x = 2exp(— ^m), thus 



e\C\~-9s-HZ/N]. 



(5) 



Remarkably, the bulk entropic contribution to the chem- 
ical potential exactly cancels the solvation energy of the 
counterion. This surprising cancellation occurs because 
X is pinned by the solvation energy gM, but will break- 
down when a significant fraction of the charged micelles 
contain counterions. Consequently, the surface potential 
depends primarily on the competition between surface 
entropy and the solvation energy of an ion on the surface. 
However, since Z depends on the ionic strength, Eq.ll^J, 
there is some residual dependence on gM- Since the left 
hand side is positive definite, meaningful results are ob- 
tained only when — ln[Z/iV] > gs- Therefore, surface en- 
tropy's contribution to the chemical potential must over- 
come the solvation energy of the surface charge. This is 
in stark contrast to thermodynamic models of charging 
in aqueous systems, where the configurational entropy of 
the surface is typically ignored tl S] . 

Combining this relation with our structural and elec- 
trokinetic measurements, we can estimate the solvation 
energy of a charge on the particle surface, gs- As- 
suming that charges on the surface are constrained to 
reside within an adsorbed micelle, then packing con- 
straints limit the number of charges to 10^. Our interac- 
tion measurements are consistent with eC^/ksT « 4 and 
Z « 100 - 1000, therefore ~ 3 - 5 fc^T, about one- 
quarter to one-half of the energy cost of charging an ion 
in bulk. 

We have found that the electrostatic interactions of col- 
loidal particles in a nonpolar solvent can be surprisingly 
strong and tunable. While the functional form of inter- 
particle potentials is identical the those of charged parti- 
cles in aqueous solvents, the thermodynamics of charging 
is qualitatively different. We expect reverse- micelle me- 
diated colloidal interactions to play an important role in 
the stability and structure of new colloidal systems and to 
enable the exploration of fundamental physical phenom- 
ena in previously unattainable regimes. Moreover, our 
framework may help to illuminate the origin and conse- 
quences of charging in technologically relevant nonpolar 
systems. 
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